A Potential Energy Landscape Study 
of the Amorphous- Amorphous Transformation in H 2 
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We study the potential energy landscape explored during a compression-decompression cycle for 
the SPC/E (extended simple point charge) model of water. During the cycle, the system changes 
from low density amorphous ice (LDA) to high density amorphous ice (HDA). After the cycle, the 
system does not return to the same region of the landscape, supporting the interesting possibility 
that more than one significantly different configuration corresponds to LDA. We find that the 
regions of the landscape explored during this transition have properties remarkably different from 
those explored in thermal equilibrium in the liquid phase. 



The physics of water in its supercooled and glass 



sv 
I, 



phases is the object of several recent experimental 
S IE 0, IE and theoretical 0, H, H E3 investigations. 
One of the most fascinating aspects is the possibility of 
a liquid-liquid phase transition 8] and the relation be- 
tween the two liquid phases and the glassy phases of wa- 
ter |9j- These glassy phases are produced via routes as 
different as crystal compression, hyper-quenching, and 
vapor deposition, and are characterized by different den- 
sities and by different local structures [j, yj [Tj. The 
intrinsic out-of-equilibrium nature of these phases leads 
to some ambiguities in the identification of the materi- 
als produced via different routes, and recent debate con- 
cerns the classification of the different amorphous struc- 
tures found [E IE El High density amorphous ice 
(HDA) resulting from the compression of hexagonal ice 
at 77 K can be converted to low density amorphous ice 
(LDA) by releasing the pressure at « 115 K [Ta . IT 3. IT ^ . 
HDA and LDA can be inter-converted via an apparent 
first-order transition by applying or releasing pressure 
[Til HEj. Molecular dynamics simulations, based on dif- 
ferent model potentials [lfjj, have been able to reproduce 
the qualitative features of the crystalline ice-HDA, HDA- 
LDA, and LDA-HDA transitions @,@> even if the com- 
pression rates in simulations are several orders of magni- 
tude larger than in experiments. 

To study the LDA-HDA transition and the relation be- 
tween glassy and liquid water we use tools that have been 
recently developed to study out-of-equilibrium liquids 

ein the potential energy landscape (PEL) framework 
In the PEL approach, the 6iV-dimensional 
configurational space — defined by the 3N center of mass 
coordinates and by the 3N Euler angles, where N is the 
number of molecules — is partitioned into a set of basins, 
each of them associated with a different local minimum of 
the potential energy landscape. The set of points belong- 
ing to the same basin are those which, under a steepest 
descent minimization procedure, end up in the same lo- 
cal minimum. The local minima in the PEL are called 
inherent structures (IS). Each IS configuration is char- 



acterized by its potential energy (e/s), its pressure (Pis) 
and by the 67V local curvatures of the PEL, which can be 
estimated in the harmonic approximation by diagonaliz- 
ing the Hessian matrix 1211. Basins with different depths 
have different curvatures |22l I23L 123. Esj . One possible 
measure of the average basin curvature — the one used 
in the present work — is offered by the "shape function" 
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where cui is the frequency of vibrational mode i \2l\ , h is 
the Planck constant, and Aq = 1 kJ/mol. This choice is 
motivated by the fact that 5is enters into the evaluation 
of the basin free energy in the harmonic approximation 

Recent work has provided a detailed statistical descrip- 
tion of the number, depth and sh ape of the basins, for 
models of both simple liquids 0, HE H3 an d molecu- 
lar liquids [25|,|28|. It has been shown that, on cooling, 
the liquid samples regions of the PEL characterized by 
lower and lower ers values |2fl and that both Pig and 
iSis are correlated with eis [2|| ■ Numerical studies of the 
PEL sampled by the equilibrium liquid permit precise 
calculations of both iSis(eis,p) and Pis(eis,p), providing 
a detailed description of the region of the PEL sampled 
under equilibrium conditions [U EE El ■ 

In this Letter we aim at comparing the properties of 
the PEL sampled during the LDA-HDA transition with 
the properties of the PEL explored by the equilibrium 
liquid. If the regions of PEL explored by the glass dur- 
ing the transition are identical [l^ to the regions ex- 
plored by the equilibrium liquid, we are entitled to con- 
nect the glass-glass transformation (GGT) to a transi- 
tion taking place in the liquid state, supporting the hy- 
pothesis that the LDA-HDA transformation is an out- 
of-equilibrium manifestation of a liquid-liquid first order 
transition. We discover that during the transformation 
the system explores regions of the landscape which are 
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never explored in equilibrium, supporting the possibility 
that the LDA-HDA transformation and the liquid-liquid 
first order transition might be independent phenomena. 
In this respect, the LDA-HDA transformation observed 
in numerical simulations 7, 8, 9] should not be taken as 
proof of the existence of a liquid-liquid transition. 

We perform molecular dynamics (MD) simulations of 
216 molecules interacting via the simple point charge ex- 
tended (SPC/E) model of water |2g. This model has 
been studied extensively; the p and T dependence of 
structural and dynamic properties in equilibrium have 
been calculated. We use a simulation time step of lfs 
and long range forces are handled using the reaction field 
method. We identify the IS by minimizing the poten- 
tial energy using a conjugate gradient minimization al- 
gorithm and calculate eis , <?is , and Pis in the resulting lo- 
cal minimum configuration. During the MD simulation, 
starting equilibrium configurations at p = 0.9 g/cm 3 are 
either (a) decompressed to 0.8 g/cm 3 or (b) compressed 
to 1.4 g/cm 3 and then decompressed to 0.8 g/cm 3 . We 
simulate at both T = K and T = 77 K. 

At 77 K we perform MD simulations using two dif- 
ferent compression/decompression rates (dp/dt = 5 x 
10 _4 g/cm 3 /ps and 5 x 10~ 5 g/cm 3 /ps) and average over 
16 different realizations. In the calculation at T = K, 
each step consists of a density change of either Ap = 
5 x 10~ 4 g/cm 3 or Ap = 5 x 10~ 5 g/cm 3 , followed by an 
energy minimization. At each step the system is com- 
pressed by Ap and the center of mass of each molecule is 
isotropically scaled. Initial configurations are extracted 
from a pre-existing ensemble of IS generated by quench- 
ing equilibrium liquid configurations at p = 0.90 g/cm 3 
and T = 220 K. At this state point, SPC/E describes the 
LDA structure accurately. Any low-temperature equilib- 
rium configuration with density close to the LDA density 
could be used as a starting configuration. 

Figure ^ shows the behavior of Pjg and eis for the 
compression and decompression of a single T = K con- 
figuration. Following each density change, the system 
is displaced from the local minimum in the PEL and a 
steepest descent minimization is performed to bring the 
system to the new location of the minimum, eis and Pis 
are continuously modified under the density changes, and 
the initial parts of the curves in Figs, andQJ) show 
such continuous modification. Above p w 1.05 g/cm 3 , 
this continuous, smooth process is interrupted by sud- 
den changes in both Pis and eis, clearly visible in Fig.^ 
Figure \Q shows the squared distance between two IS 
differing by 10~ 2 g/cm 3 , |AP| 2 = ^E^drf, where 
dri is the displacement of the center of mass of the i-th 
molecule. We see that |AP| 2 changes significantly when 
discontinuous jumps in Pig and eis occur. Thus, as in 
silica |34l |35j , these changes are due to a mechanical in- 
stability associated with the vanishing of the lowest fre- 
quency mode, which forces the system to abandon the 
previous unstable configuration in favor of a new basin. 

Figure |3 shows Pis, eis, an d <?is for both T = K 
and T = 77 K, averaged over 16 different realizations 






■f ° 

r 
I o 





003 

002 

001 

000 
0. 



Compression 
Decompression 




75 0.85 0.95 1.05 1.15 1.25 1.35 1.45 

P [g/cm ] 



FIG. 1: (a) P m and (b) e IS for T = K as functions of 
p, where a starting glass configuration (black dot) is decom- 
pressed to p = 0.8, and also compressed to p = 1.4 g/cm 3 
and then decompressed to p — 0.8g/cm 3 . The density in each 
step of the compression/decompression cycle is changed by 
Ap = 5 x 10 _5 g/cm 3 . (c) Squared distance between a se- 
quence of IS configurations differing by 10~ 2 g/cm 3 along the 
compression/decompression path. Note the correlation be- 
tween the large molecular rearrangements in part (c) and the 
corresponding discontinuities in Pig and eis in part (a) and 
(b), such as one finds at p = 1.20g/cm 3 . 



as functions of p during compression and decompression 
at two different rates. The Pis curve shows the typical 
behavior found for the pressure P observed in previous 
studies of com pres sion and decompression of tetrahedral 
glasses 0,0,0^3- We find a strong hysteresis, an obser- 
vation which has been interpreted as evidence in favor of 
a first order transition between two distinct structures, 
LDA and HDA 0,0||. 

On compressing, a significant flattening of Pis is ob- 
served when the density becomes larger than « 1.1 g/cm 3 
at T = K and larger than w 1.0 g/cm 3 at T = 77 K. At 
the beginning of the compression, the compression rate 
does not appear to play a significant role, but when the 
system is forced to change basins (due to thermal effects 
or mechanical instabilities) the compression rate becomes 
quite relevant. For the smaller compression rate, Pis and 
eis are slightly smaller. 

We next compare our results for T = K and T = 
77 K. Transitions at T = K are driven by the vanish- 
ing of the lowest normal mode frequencies resulting in 
the compression-induced disappearance of the explored 
basin. From the top panel of Fig. |3 we see that even 
though both systems start from the same configuration 
(and hence same IS) at p — 0.9 g/cm 3 , the additional 
thermal energy of the 77 K system enables the escape 
from the starting basin to occur at a smaller density. 
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FIG. 2: Pjg, eig, and <5>ig, averaged over 16 different realiza- 
tions, as functions of p during the compression and decom- 
pression at T = K (left panels) and T = 77 K (right panels) 
for the fast (solid line) and slow (dashed line) compression 
rates. Error bars are shown on the left lower corner. 



This results in a more effective relaxation process and, 
correspondingly, in a flattening of Pig at a smaller den- 
sity value (compared to the T = K case) and in smaller 
Pis and eis values. The shape of the basin is also signifi- 
cantly modified by the compression/decompression cycle. 
We also find that the basin shape £>is increases on com- 
pression and decreases on decompression, coherent with 
the frequency shifts of the translational peaks of the den- 
sity of states (i.e. in the probability distribution for u>i). 
In the region where Pis flattens, neither Sis nor the den- 
sity of states changes significantly. 

After a complete compression/decompression cycle, 
the system does not return to the same configuration. 
Indeed, Fig. shows that while Pis has the same value 
on compression and decompression at p k 0.89 g/cm 3 , 
the values for eis and <Sis differ — i.e., atpw 0.89 g/cm 3 , 
T and P are identical, but quite significant differences arc 
detected at a microscopic level. This finding supports the 
interesting possibility that more than one significantly 
different configuration corresponds to LDA. 

We next compare PEL properties sampled in equilib- 
rium at different T and p with PEL properties sampled 
by the glass during the compression/decompression path 
with the slow compression rate. Figure 3 contrasts Pis 
and Sis m & n equilibrium liquid (solid symbols) |30j and 
during the compression/decompression runs (open sym- 
bols) at T = 77 K. Surprisingly, we find that even a slight 
density change (from p = 0.9 to 1.0 g/cm 3 ) allows the 
system to begin exploring regions of configuration space 
not explored under equilibrium conditions at the corre- 
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FIG. 3: (a) Pis and (b) <Sig as functions of eis a t dif- 
ferent densities for the equilibrium liquid (filled symbols) 
and for the compression/decompression of glass (thin line) 
at both temperatures studied. Open symbols indicate the 
density during the compression/decompression. The rela- 
tive location of the open symbol with respect to the loca- 
tion of the same filled symbol is a measure of the difference 
between the landscape properties in equilibrium and in the 
glass. The significant difference in the location confirms that 
the compressed/decompressed glass explores PEL regions not 
explored under equilibrium conditions at the same density. 
Densities are 0.9(O), 1-0(0) , Ll(<), L2(v)i L3(t>), and 
1.4(A) g/cm 3 . The error bars for these densities are indi- 
cated in the figures. Data for equilibrium liquid are averaged 
over 2000 realizations and error bars are smaller than the cor- 
responding symbol size. 



sponding density. For a given p, Pis and <5>is are larger 
in the compressed glass than observed under equilibrium 
liquid conditions. We find that the T = K results show 
even larger differences. 

For small density changes (p < 1.1 g/cm 3 at T = K, 
or for p < 1.0 g/cm 3 at T = 77 K) the system is still 
confined in the (deformed) starting basin. In this den- 
sity range (i.e., before any basin change) compression rate 
effects do not play any role. Figure shows that even be- 
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fore basin change effects become relevant, the system is 
already in a compressed state which hardly remembers 
an equilibrium liquid configuration. The values of Pis 
and iSis at p = 1.0 g/cm 3 are comparable to the equilib- 
rium liquid values corresponding to a much higher den- 
sity of 1.3 g/cm 3 . Hence we conclude that the flattening 
of the Pis curve of Fig. |2 — associated with the LDA- 
HDA transition — takes place in a region of the PEL that 
is never explored under equilibrium conditions. 

Our results support the view that the compres- 
sion/decompression of glassy water takes place along a 
PEL path rarely explored by the equilibrium liquid. Such 
information — which could not have been extracted from 
standard analysis (such as standard thermodynamic or 
structural quantities) — is made possible by the ex- 
treme sensitivity of the landscape properties. We have 
shown that the IS visited by the liquid (on the time 
scale probed in computer simulations) are different than 
those sampled by the glass. This observation has pro- 
found consequences for the possibility of developing a 
thermodynamic approach to the LDA-HDA transforma- 
tion. Indeed, only when the glass state can be associated 
with a frozen liquid (i.e. when the equilibrium relations 
>5is(eis,p) and Pis (eis, p) are satisfied) it becomes possi- 
ble to provide a thermodynamic description of the glass 
and formally connect the mechanical instability of the 
LDA-HDA transformation to a mechanical instability in 



the liquid state^, |2(|. If the liquid and the glass had 
shared the same portion of the PEL, then the mechani- 
cal instability observed during the transformation would 
have constituted genuine evidence of a liquid-liquid tran- 
sition. The observation that the LDA-HDA transforma- 
tion is unrelated to properties of the equilibrium liquid 
[3fil | does not allow us to relate the observed numerical 
GGT to a liquid-liquid transition. Due to the extremely 
different time scales probed in experiments and in simu- 
lations, our numerical results do not exclude neither the 
existence of a liquid-liquid transition nor the possibility 
that the LDA-HDA transition observed experimentally 
is related to a liquid-liquid transition. Our results show 
that the PEL properties at the beginning and the end of 
the compression/decompression cycle are different, sug- 
gesting that several distinct configurations can be asso- 
ciated with LDA. 
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